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ABSTRACT 

Context, n Carinae is the colliding wind binary with the largest mass loss rate in our Galaxy and the only one in which hard X-ray 
emission has been detected. 

Aims, rj Carinae is therefore a primary candidate to search for particle acceleration by probing its gamma-ray emission. 
Methods. We used the first 21 months of Fermi/LAT data to extract gamma-ray (0.2 - 100 GeV) images, spectra and light-curves, 
then combined them with multi-wavelength observations to model the non-thermal spectral energy distribution. 
Results. A bright gamma-ray source is detected at the position of 77 Carinae. Its flux at a few 100 MeV corresponds very well to the 
extrapolation of the hard X-ray spectrum towards higher energies. The spectral energy distribution features two distinct components. 
The first one extends over the keV to GeV energy range, and features an exponential cutoff at ~ 1 GeV. It can be understood as 
inverse Compton scattering of ultraviolet photons by electrons accelerated up to y ~ 10 4 in the colliding wind region. The expected 
synchrotron emission is compatible with the existing upper limit on the non-thermal radio emission. The second component is a hard 
gamma-ray tail detected above 20 GeV. It could be explained by 7T°-decay of accelerated hadrons interacting with the dense stellar 
wind. The ratio between the fluxes of the n° and inverse Compton components is roughly as predicted by simulations of colliding 
wind binaries. This hard gamma-ray tail can only be understood if emitted close to the wind collision region. The energy transferred 
to the accelerated particles (~ 5% of the collision mechanical energy) is comparable to that of the thermal X-ray emission. 
Conclusions. We have measured the electron spectrum responsible for the keV to GeV emission and detected an evidence of hadronic 
acceleration in 77 Carinae. These observations are thus in good agreement with the colliding wind scenario suggested for 77 Carinae. 

Key words. Gamma rays: stars - X-rays: binaries - X-rays: individuals: 77 Carinae - X-rays: individuals: FGL J1045.0-5942 - 
Acceleration of Particles 



1. Introduction 

About 30 early-type stellar systems feature synchrotron radia- 
tion in the radio domain, a signature of electron acceleration (De 
|Becker|2007|). Diffusive shock ac celeration in stellar wind colli- 
sions (Benaglia & Rom ero|2003[ ), either in colliding wind bina- 
ries or OB associations, is the most likely acceleration process 



and a candidate for cosmic ray acceleration ( Axford 1981 Casse 
|& Paulp982 i. Gamma-rays, emitted by hadrons accelerated in 



stellar wind collisions, have however not yet been identified 

The gamma-ray emission expected from a colliding wind bi- 
nary increases with the stellar wind mechanical luminosity, the 
fraction of the wind enduring collision, the photon energy den- 
sity for inverse Compton emission and the matter density for 
proton-proton interaction and subsequently 7r°-decay. Located at 
a distance of 2.3 kpc (Smith 2006 ), 77 Carinae is one of the most 



massive (80-120 Mq; |Davidson & Humphreys 1997 Hillier 



et al.||200T I and brightest stellar system in the Galaxy and fea- 



tures the strongest mass loss rate known. As, in addition, its col 
liding wind region is relatively wide, 77 Carinae is a primary can 
didate for gamma-ray detection. 

The "Great Eruption" of 1843 saw 77 Carinae become the sec 
ond brightest object in the sky; it was followed by another note- 
worthy flare in 1890 (see e.g. Davidson et al.|199 9; Fernandez - 



Lajus et a l.|2010 i. The Homunculus is the name given to the ex- 



tended bipolar nebula observed around 77 Carinae: it was crafted 
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by the colossal amount of matter (10-20 M ; |Srriith et al.| 2003) 
that was ejected during the Great Eruption. The second burst led 
to the ejection of ^ 1 Mq, creating the so-called little homuncu- 
lus. 77 Carinae keeps emitting matter thanks to its powerful stellar 
winds : the mass-loss rate is believed to be 10 _4 -10~ 3 Mq yr" 1 
T9781 |Hillier et al.|[200T] |Pittard & Corcoran| 



(lAndriesse et al. 



[200^|vanBoekel et al.|2003 1. 

Even if not all questions are settled, strong evidences sug- 
gest that 77 Carinae is a binary system. For instance, radio 
( |Duncan et al.|1995 1, millimeter (mm;|Abraham et al. 2005 1, op- 
tica l (|Damineli|1996 |Damineli et al.|2000[ ), near-infrared (near- 



IR; |Whitelock et al.| f 
( |Corcoran||2005| ) observations, 



1994| |2004[ |Damineli| p96), and X-ray 
obtained over the last decades, 
have unveiled the existence of a period of ~ 2022.7 + 1.3 days 



flPamineli et al.|2008) 

The first component in the binary system is very likely a 
luminous blue variable (LBV; |Davidson & Hum phreys 1997). 
The second component is probably a late-type nitrogen-rich O 
or Wolf-Rayet (WR) star ( jlping et al.|2TjQ3j|Verner et al. 



2005) 



2001 



2007 



The semi-major axis of the orbit is 16.64 AU (Hillier etal 
and the eccentricity is very high (e ~ 0.9; Nielsen et al. 
Thus, the periastron distance r P e r i a stron is around 1.66 AU, while 
the primary star's radius R* i, although poorly constrained, is es 
timated between 0.7-1 AU (resp 
|Darriineli|19"96] ). 

The X-ray emission of 77 Carinae can be discriminated into 
two components, precisely known thanks to Chandra observa- 
tions ( |Sewardetal.|2001| ). The softer part (kT 5x ~ 0.5 keV) pre- 



Corcoran & Hamaguchi 2007 
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vails in the spectrum up to 1.5 keV. It is spatially extended, and 
is likely linked to the stellar winds colliding with the interstellar 
matter. The harder part (kTux ~ 4.7 keV) dominates over the 
2-10 keV domain. It is punctual, centered on the binary system, 
and believed to be due to the hydrodynamical shock created by 
the stellar winds of both components colliding with each other. 
The X-ray spectrum of 77 Carinae supports the colliding-wind 



binary scenar io (|Usov|1992||Stevens et al.|1992||Corco ran 2005 
|Pittard|2007| |Parkin et al.|2009) . The basic idea is that the LBV 
emits a rather slow and dense stellar wind (voo.i ^ 500 km/s, 
Mi s 2.5 x 10~ 4 M /yr; |Pittard & Corcoran|2002| ), which col- 
lides into the faster, shallower wind originating from the sec- 
ondary component (Voo,2 - 3000 km/s, M 2 - 1-0 x 10~ 5 M /yr; 
Pitta rd & Corcoran|2002| l. In this framework, the wind collision 
region represents the location of the hydrodynamical shock, and 
thus of the X-ray emitting region. 

A hard X-ray tail had first been observed towards 77 Carinae 



by BeppoSAX (Viotti et al. 



2004), and has subsequently been 



unambiguously confirmed by both INTEGRAL (Leyder et al. 
[ |2010| l and Suzaku (Sekiguchi et al. 20091 observations. 



This gives a strong evidence that the wind collision leads to non- 
thermal electron acceleration, and suggests the possibility of a 
gamma-ray detection. 

The Carina region has been observed with several gamma- 
ray experiments. Agile detected the source 1AGL J1043-5931 
(ar=161.2, £= -59.7, with an error of 0.4°) that could be related 
to 77 Carinae (Tavani et al. 2009). A 2-day y-ray flare has also 
been observed in October 2008, although its origin is unclear. 
Fermi/LAT detections where reported in the Fermi bright source 
list (BSL,|Abdo et al.|2009| ) and in the first year catalog (1FGL, 
|Abdo et al.|2010) . In the BSL, the source 0FGL J1045.6-5937 
lies at w 5.4 arcmin of 77 Carinae and is therefore not associ- 
ated with it. In the 1FGL, the source 1FGL J1045.2-5942 lies 
at 1.7 arcmin from 77 Carinae which is again too far for a for- 
mal identification. A possible association with the open cluster 
Trumpler 16, which may hold some energetic young pulsars, was 
suggested ( |Abdo et al.|2010| ). In the TeV domain, no detection 
has been reported in the vicinity of 77 Carinae. 

In the following sections, we describe our detailed anal- 
ysis of the Fermi/LAT data, and discuss the resulting multi- 
wavelength spectrum of 77 Carinae in terms of electronic and 
hadronic accelerations. 



2. Fermi/LAT observations 

The large area telescope (LAT) aboard the Fermi spacecraft is an 
electron-positron pair conversion telescope sensitive to photon 
energies from 20 MeV to > 300 GeV ( |Atwood et al.|2009| ). The 
LAT is made of 16 towers, each including a particle tracker and 
a calorimeter (allowing for the reconstruction of the direction 
and energy of the incident gamma-ray photons), and of an anti- 
coincidence shield rejecting the charged particle background. 

The data are provided by the Fermi Science Support CenteiQ 
which also delivers dedicated analysis software, called the Fermi 
ScienceTool^\ The software modules allow to select events from 
an incoming region of the sky, to compute the corresponding 
live-time exposure map and the likelihood of an input sky model, 
and to derive significance maps, source spectra and light-curves. 
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2.1. Data analysis 

We have used Fermi/LAT data accumulated in the Carina region 
during 21 months, from August 4, 2008 to April 3, 2010. The 
reconstructed gamma-ray photons were selected in a circular 
region centered on 77 Carinae with a radius of 10°, belonging 
to the diffuse clas^ and with energies between 200 MeV and 
100 GeV. This upper limit is justified by the fact that at higher 
energies, systematics of the instrument and the role of cosmic- 
ray events are not yet fully understood. The low energy thresh- 
old was chosen for several reasons. Firstly, 77 Carinae is located 
in the Galactic plane, where the galactic diffuse emission (whose 
spectrum can be described by a powerlaw with a photon index of 
2.2) is very strong. Selecting a threshold of 200 MeV decreases 
the number of background events by a factor of 4. Secondly, 
the galactic diffuse emission, whose knowledge is crucial to es- 
timate the sky model parameters at low energies, is affected by 
large systematic uncertainties, influencing the likelihood analy- 
sis. Thirdly, the size of the LAT point spread function decreases 
rapidly with increasing energy (from 100 MeV to 200 MeV, 
the 68% event containment radius decreases from 5° to 3°). As 
the analysed region encompasses a large number of sources (25 
sources are listed in the 1FGL within 10° of 77 Carinae, out of 
which 10 are within 3.5°), the ability to distinguish photons from 
different sources is important. 



2.2. Localisation 

The source 1FGL J1045.2-5942 (a=161.3053, <5=-59.7057) was 
not formally associated with 77 Carinae since the latter lies 
slightly outside of the 95% confidence region. As we are using 
almost twice more data than was available when the 1FGL was 
generated, a more accurate source localisation can be obtained. 
After an initial fit of the parameters of the sources close to the 
location of 77 Carinae, as well as of the normalisations of the dif- 
fuse galactic and isotropic emissions, we determined the position 
of the nearest source to 77 Carinae with the tool gtf indsrc and 
found 0-= 16 1.265 and £=-59.7015, with a 95% confidence ra- 
dius of 1.18 arcmin. 77 Carinae is thus perfectly compatible with 
our improved position, as illustrated in Fig. [T] 

Since the width of the point spread function for photons con- 
verted in the front part of the LAT is smaller than for these con- 
verted in the back, we also derived the position using the front 
events only, but the uncertainty did not decrease significantly. 
We also derived the source location for low (< 8 GeV) and high 
(>10 GeV) energy events and obtained error circles also com- 
patible with the position of 77 Carinae (see Fig. [TJ. 

The presence of the source is indubitable, with a test statistic 
(TS, |Wilks|1938] ) > 2800 (~ 53 <x) for the 200 MeV to 100 GeV 
energy range. The low and high-energy components have respec- 
tively TS of 2281 (~ 47 cr) and 73 (~ 8.5 &). In the rest of the 
paper, we will name the Fermi/LAT source FGL J 1045. 0-5942, 
owing to its improved position. 



2.3. Spectral analysis 

The spectral analysis of FGL J 1045.0-5942 was performed us- 
ing the maximum likelihood method. The region modeling in- 
cludes two diffuse emissions (Galactic plane and isotropic) and 
35 point-like sources listed in the 1FGL catalogue. While most 
of these sources have been modeled as pure power-laws, four 
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Fig. 2. Top: Photon spectrum (left) and spectral energy distribu- 
tion (right) of FGL J1045.0-5942 with 21 months of Fermi/LAT 
data. Bottom: Spectral energy distribution for the the complete 
21 month period (red) and for the periastron period (blue). 
Uncertainties are 95% confidence levels. 



Fig. 1. Top: Fermi/LAT smoothed count map between 200 MeV 
and 100 GeV in galactic coordinates. The bright source in the 
centre of the field is FGL J1045.0-5942. The source located 
~ 1.3° to the North is the pulsar PSR J1048-5832. Bottom: 
Fermi/LAT significance map corresponding to a zoom on the 
black box shown in the upper image. The blue cross indicates 
the optical location of r\ Carinae. The green circle is the 95% 
confidence region reported for 1FGL J1045. 2-5942 (based on 
11 months of data). The plain black circle shows the 95% con- 
fidence region found in our analysis of the 21 month dataset. 
The dashed circles show the 200 MeV to 8 GeV and 10 GeV 
to 100 GeV confidence regions. The cyan circle corresponds to 
the confidence region of the hard X-ray detection obtained with 
INTEGRAL ( |Leyder et al.|20T0] i. 



among them are known to be pulsars and have thus been mod- 
eled as power-laws with exponential cutoff (Abdo et al. 2010). 
This is particularly important for the pulsar PSR J1048-5832, 
located ~ 1.3° away from FGL J1045.0-5942. 

As the detection of FGL J1045. 0-5942 is very significant, a 
spectrum could be extracted in 9 spectral bins. In each bin, the 
TS is larger than 25 (5 cr). The resulting photon spectrum and 
spectral energy distribution are presented in Fig. [2] Both binned 
( |Cash|1979] l and un-binned ( |Mattox et al.|1996| > likelihoods have 
been used and their results are in full agreement. A curvature 
at low energy and a modification of the spectral slope at high- 
energy are clearly detected. 

The overall spectrum of FGL J1045.0-5942 differs from the 
standard power-law used to build the 1FGL catalogue, and has 
been modeled as a combination of two components: a power- 
law of photon index F with an exponential cutoff plus another 



power-law at the highest energies. Since only a limited num- 
ber of spectral shapes are available in the Fermi ScienceTools, 
the dataset has been split into two energy ranges (0.2-8 GeV 
and 10-100 GeV) to determine the two spectral components in- 
dependently. We maximized the sum of the likelihood in both 
spectral bands together so that the parameters of all sources are 
constrained by the complete dataset. The resulting parameters 
are T = 1.69 ± 0.12, £ cut = 1.8 ± 0.5 GeV, F . 2 -iooGev = 
1.52 x 10~ 7 ph cirT 2 s _1 for the exponentially cutoff power-law 
and T = 1.85 ± 0.25, F . 2 -ioo GeV = 0.41 x 10~ 7 ph cirT 2 s" 1 for 
the high-energy component. 



2.4. Timing analysis 

We constructed the light-curve of FGL J 1045. 0-5942 using the 
maximum likelihood technique with time bins of one month (see 
Fig. 13) . The light-curve is consistent with a steady flux. The 
time-averaged flux of the source in the 0.2-100 GeV energy 
range is Fo.2-100 GeV = (1-93 ± 0.05) x 10~ 7 ph cirT 2 s -1 . 

The flux was also determined for the period when rj Carinae 
was at the X-ray minimum, as measured by RXTE, (i.e. from 
MJD = 54843 to MJD = 54883). The flux level obtained for 
this period is reported in red in Fig. [3] It is compatible with the 
average flux at the 95% confidence level. The spectrum detected 
during the periastron period is also compatible with the average 
spectrum obtained over the 21 months period (bottom of Fig. [2}. 
In particular the cutoff energy did not change significantly, nor 
the high-energy component, clearly detected above 40 GeV. 

The gamma-ray emission does not vary by more than 50%, 
much less than observed in the X-ray band at periastron, where 
a variability factor larger than 10 is observed. 
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Fig. 3. Top: Monthly light-curve of FGL J1045.0-5942 derived 
from the 21 month Fermi/LAT dataset. The red point corresponds 
to the flux level observed during the X-ray minimum. The black 
line represent the best fit result (with its error) obtained assuming 
a steady flux for the source. Uncertainties are 95% confidence 
levels. Bottom: Residuals obtained from the steady flux fit. 
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Fig. 4. The black line shows the RXTE X-ray light-curve of 
77 Carinae with the minimum close to the 2009 periastron pas- 
sage. The red points correspond to the high-energy events (E > 
20 GeV) detected by Fermi/LAT within 1° of 77 Carinae (note the 
event energy scale on the right). 



Fig. [4] shows the high-energy (> 20 GeV) events superim- 
posed on the RXTE X-ray light-curve^] Despite the fact that this 
high-energy component is very significant (TS > 110, <x > 10), 
there is also no apparent correlation between the high-energy 
component and the orbital phase of 77 Carinae. 

We also searched for gamma-ray flares using a flux aperture 
photometry method, but no significant deviation has been found. 
In particular, the two-day period during which Agile reported a 
flare from the Carina region has been investigated, but no sign 
of activity could be detected. During these two days, the average 
flux detected by Fermi/LAT 'is F - (1.9 + 0.6) x 10~ 7 ph cm" 2 s _1 , 
well below the flux detected by A sile which is (27.0 ± 6.5) x 



10 7 ph cm s 1 (Tavani et al. 2009 1. As Fermi and Agile have 



different orbits and as we used an energy threshold of 200 MeV, 
our analysis cannot exclude a very short and low energy flare. 



3. Discussion 

3.1. Association with 77 Carinae 

Since the spectral energy distribution observed by Fermi/LAT is 
made of two components (Fig. [2]), a spatial coincidence of two 
distinct gamma-ray sources (for instance a pulsar and a blazar) 
cannot be excluded a priori. However this coincidence should 
match with the centroids of the low- and high-energy compo- 
nents, which are separated by less than 1.2' (Fig.[T| and compat- 
ible within the uncertainties. This is rather unlikely as this would 
require two bright gamma-ray sources very close to 77 Carinae. 

Given the number of pulsars detected within 2° of the 
Galactic plane (26 or 28^]), the probability to find one of them in 
the 95% confidence region of 77 Carinae is 2.5 x 10~ 5 . Moreover, 
we conducted a blind search frequency analysis and did not find 
any evidence for a pulsation in FGL J1045. 0-5942. 

Given the 685 sources associated with blazars in the 1FGL 
over the entire sky, the probability to find a blazar at the location 
of 77 Carinae is 2.3xl0~ 5 . In addition, if the high-energy emission 
was related to a blazar, flaring episodes would be expected, but 
none is observed (see Fig.|4|. 



4 ftp : //legacy . gsf c . nasa . gov/xte/data/archive/ 
ASMProducts/definitive_ldwell/lightcurves/ 

3 The 1FGL quotes two potential y-ray millisecond pulsars within 2° 
of the Galactic plane. 



These two results imply that the probability to find a pulsar 
and a blazar at the location of 77 Carinae is lower than 10~ 9 . 

With the above probabilities and the absence of any other 
bright X-ray source in the INTEGRAL error circle (Leyder et al. 
2010), we can safely assume that the high-energy emission de- 
tected by Fermi and INTEGRAL comes from a single source, 
very likely 77 Carinae. 



Chandra images of 77 Carinae ( Weis et aL| 2004 1 show an 
extended soft shell-like component (corresponding to an exter- 
nal shock at the boundary of the Homunculus nebula) around 
the harder, point-like emission. |Ohm et"aL ( 2010[ ) have argued 
that the gamma-ray emission detected by Fermi could be emitted 
by this external shock. Their conclusion was based on the first 
eleven month spectrum, which did not exhibit the two spectral 
components presented above. The high-energy spectral compo- 
nent cannot be explained by the external shock, because the den- 
sity is far too low to allow significant hadron interactions and n° 
cooling. In addition, |Leyder et al.| ( |2010| l have shown that the 
contribution of this outer shell to the hard X-ray emission from 
77 Carinae was lower than 15%. Thus the outer shell interpreta- 
tion seems unlikely. 

3.2. Particle acceleration in 77 Carinae 

We will assume below that particle acceleration takes place in 
77 Carinae through diffusive shock acceleration in the colliding 
wind region of the binary system. In the shock region, particle 
acceleration is counterbalanced by four main cooling processes: 

1 . Inverse Compton scattering of electrons in the intense ul- 
traviolet radiation field of the stars. The cooling time scale 
is 



tic 



3ym e c 



3nR 2 m e c 2 



4o- T cy 2 B 2 U ra d o- T jB 2 L Tio^-io'l, 



R 2 

ll)14cm x6 10 2 s, 



where L is the stellar ultraviolet luminosity of the primary 
star and R is its distance to the colliding wind region. 

2. Proton-proton interactions and subsequent pion decay. 
|Bednarek| ( |2005[ ) pointed out that high-energy hadrons in 
the wind of massive stars will be photo-disintegrated into 
protons and neutrons on a short timescale. The interaction 
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timescale for protons is inversely proportional to the density 
of matter in the post-shock region, where the particles are 
trapped by the magnetic field. Hydrodynamic simulations 
( Pittard 2009 ) indicate that the density in the shock region 
is much larger than the unperturbed wind density, by a fac- 
tor of 5 ~ 1-100. The p-p cooling time scale can therefore be 
written : 



1 

a-p p 5nc 



AnR) h m p V w 
<T pp dMc 



<5io M w - 



x4- 10 5 s, 



4 M /yr 



where V is the typical wind velocity and M is the mass loss 
rate of the primary star. 

Proton diffusion away from the shock region. The diffu- 
sion timescale is related to the bulk velocity of the post- 
shock material from the central regions towards the outside. 
Hydrodynamic simulations indicate that the material flows 
out of the shock region at a velocity smaller than, but of 
the order of, the pre-shock wind velocity (Pittard 2009}. As 
the shock region has a size similar to the stellar separation, 
the bulk diffusion timescale can be estimated as ( |Bednarek| 
[2005] l : 

t bu i k = 3R/Vr- Rw 



x3 



10 6 s. 



Vl0 3 km/s 

4. Electron bremsstrahlung in the post-shock region. For a 



density 



8n — 6 



M 



— — x 3 ■ 10 cm , 



4nR 2 Vm p Rw^cm 

the cooling timescale is given by ( |Aharonianj2 004 ) : 

tbr ~ 1.2 x 10 15 ((J« cm -0 _1 s ~ 4 • 10 s s. 

Particle acceleration in the shock is therefore mainly coun- 
terbalanced by Inverse Compton scattering for electrons and 
by proton-proton interactions. Equating the acceleration time 

face = T'(f) respectively with tjc and t pp for electrons and 
protons (Ri is the Larmor radius) provides the maximum char- 
acteristic energy of the particle distribution. 

For electrons, a powerlaw spectrum is expected, with an ex- 
ponential cutoff at an energy 



Inec 2 BR 2 V 



IG ' "l0 14 cm 



J 5-10 6 L„ 



Vl03kn,/sX3-10 4 , 



where B is the magnetic field in the shock region. 

Jmax.e should be fairly independent of the orbital position in 
the dipole approximation for a magnetic field varying as R 2 . 
The gamma-ray spectrum will therefore feature an exponential 



cutoff at the maximal energy as derived by |Eichler & Usov 
( fT993"] >. A cutoff energy of ~ 1 GeV (Sect. |23j ) corresponds to a 
magnetic field of ~ 50 G at the stellar surface. 

For protons, the maximum characteristic energy is limited by 
proton-proton interactions to 
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Fig. 5. Spectral energy distribution of 77 Carinae includ- 
ing BeppoSAX/MECS (blue), INTEGRAL/ISGRI (red) and 
Fermi/LAT (purple) data, and a radio upper limit to the syn- 
chrotron emission (green). From low to high energies are shown 
the synchrotron, stellar emission, inverse Compton and jT -decay 
spectral components. 

Table 1. Model parameters of the shock region determining the 
non-thermal spectral components of 77 Carinae, as displayed in 
Fig.0 





Parameter 


Value 




Environment 


Photon energy density 
Magnetic field 
Density 


2.7 
0.5 
3 ■ 10 9 


erg/cm 3 
G 

cm -3 


Electron distribution 


Powerlaw index 

Tmax.e 

Total energy 


2.25 
10 4 
10 40 


erg 


Proton distribution 


Powerlaw index 

ymax.p 

Total energy 


2.25 
10 4 
1.3 • 10 40 


erg 



parameters of the model are listed in Table [T] The magnetic 
field and the electron energy distribution were adjusted to match 
the upper limit on the radio synchrotron emission, derived from 
the minimal thermal emission detected with ACTA ( |Duncan &| 
White 2003 ), and to match the inverse Compton continuum de- 
termined by INTEGRAL and Fermi. The slope and cutoff en- 
ergy of the interacting proton energy distribution were fixed to 
2.25 and 10 4 , respectively, and its normalization was fitted to 
match the high-energy gamma-ray tail using the 7r°-decay model 
of Kelner et al. (2006). The proton and ultraviolet photon energy 
densities in the shock region were fixed to the values expected 
for an average distance R = 10 14 cm. 



AnR 2 eB IV 



rmax,p 



PP' 



6M 



(7)" 



R W"cm B ^ 



V 3 

10 3 km/s 



£lO-MlO- 4 M /yr 



X 4 • 10 4 



The non-thermal spectral energy distribution of 77 Carinae 
has been represented with a model consisting of two cutoff pow- 
erlaw distributions for the electrons and for the interacting pro- 
tons. The model and the data are shown together in Fig. [5] The 



3.3. Energetic Considerations 

The wind momentum ratio of 77 Carinae is 77 = (M2V2)/(Mi Vi) ~ 
0.2 ([Pittard & Corcoran"||2002| ). The half-opening angle of the 
shock region is therefore ~ 1 rad, the fraction of the wind in- 
volved in the wind collision region is ~ 10%, and the mechanical 
energy available in that region is ~ 200 L . 
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For a density of cold protons of 3 • 10 9 cm 3 in the shock with lower stellar winds or magnetic fields will feature fainter 



region, the normalization of the /r°-decay spectrum (Fig. [5]) re 
quires a total interacting proton energy E p ~ 10 40 erg. The en 
ergy injected in the shock to sustain the observed proton distri- 
bution is of the order of E p /t pp ~ 10 L Q . 

The gamma-ray spectrum of 77 Carinae thus indicates that 
~ 5% of the shock mechanical energy (or less than 1% of the 
total wind mechanical luminosity) is transferred to accelerated 
protons downstream. This is in agreement with recent numer- 
ical simulations of relativistic collisionless shocks (Spitkovsky 
2008| l, if a significant fraction of the hadrons interacts and gen- 
erates 7T°. 

Modeling the thermal X-ray emission depends on detailed 
hydrodynamical simulations (see e.g. Pitt ard & Corcoran[2 002 ) 
and could be affected by many phenomena (Par kin et al.|2 009). 
Such simulations can explain an X-ray luminosity of « 20 L Q 
above a few keV. Our analysis hence indicates that the fraction of 
the shock energy accelerating protons is similar to that emitting 
observable X-rays in 77 Carinae. 

The ratio between the inverse Compton or 7r°-decay and the 
X-ray emission observed in 77 Carinae is larger than predicted by 
existing models applied to WR 140 (Pittard & Dougherty 2006 



Reimer et al. 2006). The high-energy emissivity predicted for 
WR140 is however not well constrained and varies by a factor 
of 100 depending on the model parameters. The strong inverse 
Compton emission and the enhanced 7r°-emission (when com- 
pared to the X-ray emission) might be related to the strong ul- 
traviolet photon field and to the very high density in the wind 
collision region of 77 Carinae, strengthening simultaneously the 
proton-proton collision rate and the absorption of the thermal 
emission. Detailed modeling is outside the scope of this paper, 
but these observations of 77 Carinae provide the first observa- 
tional constraint on the fraction of the mechanical power injected 
into particle acceleration. 

3.4. On the detection of neutrinos 

The high-energy gamma-rays observed in 77 Carinae are likely 
produced by 7r°-decay, a process which produces as many neu- 
trinos as gamma-rays. Thus, the detection of neutrinos would 
provide a conclusive evidence that hadronic acceleration is at 
work in 77 Carinae. 

The neutrino spectrum resulting from the decay of n° has 
been studied by Kappes et al. (2007 1. In the case of 77 Carinae, 
and based on y ma x,p~W > ', the expected spectrum is 



dNy ^ I Ey V 



-2.35 



- V£,70.25 TeV 



x 10~ 12 (cm 2 s TeV)" 



The atmospheric neutrino background is far too high to ex- 
pect any detection at such a low energy. However, since we 
have not measured the high-energy cutoff of the 7r°-decay emis- 
sion, we cannot exclude that Cherenkov experiments might de- 
tect 77 Carinae at energies higher (y m ax,p Z 10 5 ) than anticipated. 
Unfortunately, even in that case, only a few neutrinos could be 
detected by KM3Net in 5 years of observations, not enough for 
a source detection. 

3.5. Other sources 

The density and magnetic field in the vicinity of the shock in 
77 Carinae are such that a large fraction of the energy carried 
by high-energy protons is converted to gamma-rays. Systems 



n -decay emission. 

Scaling the emissivity of 77 Carinae by the mass loss rate (i.e. 
the wind density) suggests that LBV systems could be detected 
by Fermi up to the Galactic center while WR or OB systems 
could be detected within 1 kpc and 0. 1 kpc, respectively. Very 
few objects could therefore be detected. OB associations could 
have a combined stellar mass loss rate as strong as 77 C arinae 
and similar efficiency in creating n° (Ozernoy et al. 1997| . They 
could thus be detected up to distances of several kpc, although 
their extended emission is more difficult to measure, especially 
in the Galactic plane. 

The mechanical energy in the stellar wind of a massive star, 
integrated over its lifetime, could reach some 10 50 erg. From an 
energetic point of view, extrapolating on the particle acceleration 
efficiency of 77 Carinae, stellar winds ejected during the massive 
star evolution in a suitable environment could be as effective as a 
SNR to accelerate hadrons. Future Cherenkov observations will 
tell if stellar wind collisions could accelerate particles up to to 
the knee of the cosmic-ray spectrum. 

4. Conclusions 

We have detected a bright gamma-ray source at the position of 
77 Carinae using 21 months of Fermi/LAT data. Its flux at a few 
100 MeV corresponds very well to the extrapolation of the hard 
X-ray spectrum of 77 Carinae (as measured by INTEGRAL and 
Suzaku) towards higher energies. The spectral energy distribu- 
tion, which corresponds to an average over almost half of the 
orbit of 77 Carinae, features two spectral components. 

The first one is a powerlaw extending from keV to GeV en- 
ergies, with an exponential cutoff at ~ 1 GeV. This component 
can be understood assuming inverse Compton scattering of stel- 
lar photons by electrons accelerated up to y ~ 10 4 in the wind 
collision region. The observed cutoff energy implies a magnetic 
field ~ 50 G at the stellar surface. 

The second component, a hard gamma-ray tail, is detected 
above 20 GeV. This bright component could be explained by 
7T° -decay of accelerated hadrons interacting with the dense stel- 
lar wind in the shock region. The ratio between the fluxes of 
the 7T° and inverse Compton components is roughly as predicted 
by simulations (Pittard & Dougherty 2006; Rei mer et al.|2 006). 
Bremsstrahlung emission is expected at a much lower level and 
with a cutoff energy similar to that of the inverse Compton com- 
ponent. 

The hard gamma-ray tail can only be understood if emitted 
close to the wind collision region. Indeed, the external shock 
between the Homunculus nebula and the interstellar medium 
occurs at a density and magnetic field strength that are much 
too small to explain the observed emission. It is possible that 
a part of the soft gamma-ray emission comes also from inverse 
Compton scattering of infrared photons in the external shock as 
suggested by |Ohm et al. ( 2010| > (see also |Leyder et al.|2010 1. It 
is however difficult to account for the hard tail without a signif- 
icant fraction of the soft gamma-ray emission being emitted by 
the colliding wind region. Detection of any hard X-ray or soft 
gamma-ray variability would rule out the external shock model. 

The energy transmitted to the accelerated particles (~ 5% of 
the wind collision mechanical energy) is of the same order as 
that of the observed thermal X-ray emission, providing an im- 
portant observational constraint for future numerical hydrody- 
namical models of the colliding wind region in 77 Carinae. 

With the exception of the short flare detected by Agile, the 
averaged hard X-ray (Leyder et al. 2010) and soft gamma-ray 
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emissions do not vary by more than 50% along the orbit, not 
even during the periastron passage. This provides an interesting 
constraint on the evolution of the magnetic field and of the wind 
collision region along the orbit. 

We finally suggest that Cherenkov telescopes could measure 
the hadronic cutoff energy in 77 Carinae. 
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